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ABSTRACT 

We present Far Ultraviolet Spectroscopic Explorer (FUSE) observations of the 
young, compact planetary nebula (PN) SwSt 1 along the line of sight to its central 
star HD 167362. We detect circumstellar absorption lines from several species 
against the continuum of the central star. The physical parameters of the nebula 
derived from the FUSE data differ significantly from those found from emission 
lines. We derive an electron density n e = 88OOI2400 cm_3 from the column density 
ratio of the excited S III fine structure levels, which is at least a factor of 3 lower 
than all prior estimates. The gaseous iron abundance derived from the UV lines 
is quite high ([Fe/S] = —0.35 ± 0.12), which implies that iron is not significantly 
depleted into dust. In contrast, optical and near- infrared emission lines indicate 
that Fe is more strongly depleted: [Fe/H]= -1.64±0.24 and [Fe/S]= -1.15±0.33. 
We do not detect nebular H 2 absorption, to a limit iV(H 2 ) < 7 x 10 14 cm -2 , at 
least four orders of magnitude lower than the column density estimated from 
infrared H2 emission lines. Taken together, the lack of H2 absorption, low n e , 
and high gaseous Fe abundance derived from the FUSE spectrum provide strong 
evidence that dense structures (which can shield molecules and dust from the 
destructive effects of energetic stellar photons) are not present along the line 
of sight to the central star. On the other hand, there is substantial evidence 
for dust, molecular material, and dense gas elsewhere in SwSt 1. Therefore, we 
conclude that the nebula must have an inhomogeneous structure. 

We detect nebular absorption at 1040.94 and 1041.69 A from the two ex- 
cited fine structure levels of neutral oxygen. These levels give rise to far-infrared 
emission lines at 63 and 145 /im which are often used to infer gas properties, par- 
ticularly temperature, under the assumption that they are collisionally excited. 
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We find that the O I fine structure levels in SwSt 1 have an inverted popula- 
tion ratio. This requires a non-thermal excitation mechanism, which we identify 
as fluorescent excitation by the stellar continuum. To the extent that fluores- 
cence affects the level populations, the far-infrared [O I] line strengths cannot be 
directly used as diagnostics of density and temperature. 

Subject headings: planetary nebulae: general — planetary nebulae: individual 
(SwSt 1) — stars: AGB and post-AGB — nucleosynthesis, abundances — ultraviolet: 
ISM — dust, extinction 



1. INTRODUCTION 

The planetary nebula (PN) SwSt 1 (Swings & Struve 1940) is believed to be in the 
incipience of this stage of evolution, because of its small angular diameter (Kwok et al 1981; 
Aaquist & Kwok 1990) and high electron density, found by various studies to have a value 
in the range n e ~ 3xl0 4 to 10 5 cm -3 (Flower et al. 1984; de Freitas Pacheco & Veliz 1987, 
hereafter dFP87; De Marco et al. 2001, hereafter DM01). The low excitation of its nebular 
emission line spectrum and cool temperature of the central star, T e s = 40,000 K (DM01) 
are also suggestive of youth. The central star of SwSt 1, HD 167362, has been classified as 
spectral type [WC9] (Crowther et al. 1998), a class of H-deficient, C-rich PN central stars 
that exhibit emission lines indicative of strong winds, although its emission lines are weaker 
than typical members of this class (Tylenda et al. 1993; Gorny & Tylenda 2000). While the 
evolutionary origin of this class of stars is still unclear, it may be related to late thermal 
pulses (e.g. Herwig et al. 1999; Blocker 2001) or binary interactions (De Marco & Soker 
2002). 

In this paper, we present measurements of nebular absorption lines against the con- 
tinuum of the central star of SwSt 1. The Far Ultraviolet Spectroscopic Explorer (FUSE) 
spectrum probes material along a narrow "pencil-beam" toward the central star, and there- 
fore samples a very small portion of the nebula compared to emission line studies. The 
FUSE waveband includes the resonance Lyman and Werner bands of molecular hydrogen, 
several lines of Fe II and Fe III (which we use to derive the Fe abundance in SwSt 1), and 
Ge III A1088.46, which has been used previously to derive the Ge abundance of five PNe 



1 Bascd on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE 
is operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985. 
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(Sterling et al. 2002, hereafter SDB; Sterling & Dinerstein 2003). This spectral region also 
contains several other lines of species in the H II and neutral atomic regions. 

The structure of this paper is as follows: In §2, we describe the observations and data 
reduction. We discuss the method of analysis in §3, and present the Voigt component fits 
used to derive the column densities, central velocities, and Doppler spread parameters of 
several nebular species. In §4, we use these results to investigate properties of the ionized 
gas, such as the electron density and gas phase iron abundance. We also examine neutral 
and molecular gas in SwSt 1 (§5), deriving an upper limit to the column density of nebular 
molecular hydrogen in the line of sight, and investigating the excitation of the O I fine 
structure levels. The implications of our study are discussed in §6, and the results are 
summarized in §7. 

2. OBSERVATIONS 

The FUSE satellite provides access to the wavelength range 905-1187 A at a nominal 
resolution R = 15,000 (FWHM = 20 km s" 1 ). The design (Moos et al. 2000; Sahnow 
et al. 2000) consists of four prime focus telescopes and Rowland-circle spectrographs with 
microchannel plate detectors. Two of the channels are coated with Al:LiF, which provides 
optimal throughput for the spectral range 1000-1187 A, while the other two have SiC coatings 
for optimal reflectivity below 1000 A. Each detector consists of two segments, so that eight 
channel segments provide coverage of the entire spectral range. Spectral redundancy ensures 
that nearly all wavelengths are covered by multiple segments. 

SwSt 1 was observed in 2001 August with the large (LWRS, 30"0 x 30"0) aperture in 
time-tag mode for a total of 39,587 seconds. These data are available at the Multi-Mission 
Archive at the Space Telescope Science Institute (MAST) under the data set identifications 
B0690101 and B0690102. We used CALFUSE (v2.0.5), the standard reduction pipeline 
available at Johns Hopkins University, to calibrate the raw data. The exposures were then 
cross-correlated and co-added in order to increase the signal-to-noise ratio (S/N). We re- 
binned the data by three pixels, which improved the average S/N to ~15-25 per rebinned 
sample, depending on the region of the spectrum. This increased the effective pixel size from 
~ 7 mA to ~ 20 mA, corresponding to a velocity sampling interval of ~ 6.0 km s _1 near 
1000 A and ~ 5.2 km s" 1 near 1150 A. Figure 1 shows an excerpt from the FUSE LiF 2A 
spectrum of SwSt 1 from 1120-1150 A. 

We checked for terrestrial airglow contamination of the spectrum (which can be impor- 
tant for lines such as O I) by comparing the exposures taken during orbital night to the 
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full dataset. The night spectrum is essentially indistinguishable from that including daytime 
data, and we therefore conclude that airglow emission has a negligible effect on our results. 

We also employ archival Hubble Space Telescope (HST) ultraviolet observations of 
SwSt 1, acquired by the Space Telescope Imaging Spectrograph (STIS; Woodgate et al. 
1998; Kimble et al. 1998). This data was obtained in 2000 September, as part of Guest 
Observer program GO8304 (principal investigator S. Kwok). SwSt 1 was observed for a 
total of 2300 seconds with the Qf!2 x 0'.'2 aperture and E140M echelle grating, which provides 
a spectral resolution of R = 45,800 (FWHM ~ 6.5 km s _1 ) in the spectral region 1 150 
1730 A. The data was reduced with the CALSTIS v2.11 processing pipeline (Lindler 1999). 
Unfortunately, the E140M spectrum has a low S/N, particularly below 1200 A (S/N < 8), 
and therefore we make limited use of this data. 



3. ANALYSIS 

3.1. Line-Fitting Procedure 

We analyze the FUSE data in the segments with the best S/N (i.e. largest effective 
area) for a given wavelength region: LiF 1A for 1000-1083 A, and LiF 2A for A > 1088 A. 
Redundant coverage in the other spectral segments is used to confirm the LiF 1A/2A results. 
We normalize each absorption profile by fitting a low-order Legendre polynomial to the con- 
tinuum and dividing the spectrum by this function. Atomic wavelengths and data are taken 
from Morton (2000; 2003), while H_2 wavelengths and oscillator strengths are from the compi- 
lation of lines by McCandliss (2003; see also http:/ /www. pha.jhu.edu/~stephan/H2ools/hlh2data/, 
based on Abgrall et al. 1993a,b). 

The line of sight toward the central star of SwSt 1 passes through two interstellar clouds, 
as well as the absorbing nebular material. In order to separate the different absorption com- 
ponents of lines, we perform Voigt component fits to each line profile. These fits yield 
detailed information about the nebular features, such as central velocity < v >, Doppler 
spread parameter b, and column density N. The model profiles are convolved with an in- 
strumental spread function, a Gaussian of FWHM = 20 km s" 1 for FUSE, and 6.5 km s -1 for 
the STIS E140M data, to account for the systematic broadening effect from the instrumental 
resolution. 

In many cases, more than one line from a given species can be measured. When possible, 
we perform simultaneous fits to absorption lines originating in the same level, in order to 
better constrain the fits and reduce systematic errors. In a simultaneous fit, the parameters 
< v >, b, and N are forced to be the same for all lines from a given level (although the 
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central velocities may differ by a small offset, particularly for lines measured by different 
instruments). Simultaneous fitting is especially important for the medium- resolution FUSE 
spectrum, where saturation may not be evident due to instrumental broadening. We use 
this technique in several cases to test the importance of saturation effects, and to constrain 
the parameters of mildly saturated lines which, if modeled individually, would otherwise 
be highly uncertain. This procedure is most useful when the fitted lines have significantly 
different /-values, and at least one is not saturated. 

We use a different line fitting method in this paper than in previous work for the same 
dataset (SDB). While we measure several more nebular lines in SwSt 1 than SDB reported, 
our two studies have lines of S II, S III and Ge III in common. By comparing our Table 1 with 
Table 1 of SDB, one may note that the column densities we measure here are systematically 
higher than those of SDB, by about 0.3 dex (a factor of two). This is not due to differences in 
the atomic data used; indeed, the /-values for S III AA1015.78, 1021.32 are larger in Morton 
(2003) than Morton (1991). The difference is likely due to the fact that we account for 
saturation effects in the present study. We consider the results presented in this paper to 
be more robust than those of SDB, since we utilize simultaneous fitting, and more recent 
atomic data (Morton 2003) and solar abundances (Asplund et al. 2004). 



3.2. Line Fit Results 

Given the low excitation of the SwSt 1 emission line spectrum (DM01; dFP87), it is not 
surprising that the only nebular species we detect in the FUSE spectrum are neutrals and 
low ions: Ar I, O I, Fe II, Fe III, Ge III, and S III. Of these ions, S ++ requires the highest 
energy to create (23.3 eV), which is quite modest by PN standards. There is no evidence of 
absorption from more highly ionized species. 

In Table 1, we show the results of the line fits for all detected nebular species. Two lines 
are listed in Column 1 when a simultaneous fit was performed, while one line is recorded 
otherwise. The quoted l-a error estimates include continuum placement and statistical 
fit uncertainties. We estimated the statistical errors with Monte Carlo simulations of the 
absorption line profiles - 100 artificial spectra were generated for each profile, with the data 
points modified within the errors of the Gaussian noise. In most cases, the systematic errors 
from continuum placement dominate the uncertainties. 

For lines arising from different levels of the same ion, we further constrain their fits 
by forcing the lines to have the same < v > and b, although the column densities of the 
levels are allowed to differ. This is a useful constraint for ions which exhibit absorption lines 
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from several excited levels (Fe III, O I, and S III), some of which may be saturated. Since 
the excitation energies of the ground term levels for these three ions (< 0.13 eV) are small 
compared to the gas temperature (0.9 eV for T = 10,500 K; DM01), the same value of b 
should apply for lines from all the observed levels of these ions. 

The results of Table 1 are presented graphically in Figures 2a and 2b, where the 
continuum-normalized profiles are displayed versus heliocentric velocity (all velocities re- 
ported in this paper are in the heliocentric frame). Note that the line of sight toward SwSt 1 
intersects two interstellar components, one near —5 km s _1 , and the other at +20 km s" 1 
(these are seen most clearly in the profile fits to Ar I and Fe II) . These absorbing clouds are 
well-separated in velocity from nebular absorption. The profiles of S II AA1250.58, 1253.81 
and O I AA1304.86, 1306.83 are from STIS data, while all others are from the FUSE spec- 
trum. The total model fit to each line profile (convolved with the instrumental spread 
function) is shown as a thick, solid line; thin, dashed lines represent individual component 
fits before instrumental broadening. The depth of the dashed lines therefore indicate the 
importance of saturation for each component. In Figures 2a and 2b, the panels are grouped 
by species - double panels for simultaneous fits, and single ones for single line fits. The 
absorbing ion is indicated above the panel 2 , and the wavelength of the line is shown within 
the panel. Nearby absorption from other species is indicated in the figures. 

The central velocities we measure agree well with the findings of Dinerstein et al. 
(1995b), who used high resolution observations of the Na I AA5889, 5895 doublet to find 
a heliocentric velocity of v Q = —47 km s _1 for the circumstellar gas. This agreement is 
somewhat surprising, since ionized species will not necessarily absorb at the same velocity 
as neutral sodium, which has an ionization potential of only 5.1 eV. Note that the species 
with the most aberrant nebular velocities, Ar I and O I, are blended with a heavily saturated 
interstellar component, and < v > is highly uncertain. 



4. RESULTS FOR THE IONIZED REGION 

4.1. Sulfur Ionic Column Densities 

Sulfur absorption lines are useful for calibrating the abundances of other elements, as 
well as for determining the average electron density n e along the line of sight. While the 
excited levels of the S III 3 P ground term are scantly occupied in the diffuse interstellar 



2 Asterisks denote absorption from excited levels; e.g. Fe III** is from the second excited level above the 
ground state in Fe III. 
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medium (ISM), their populations can be significant in the physical conditions of PNe. Neb- 
ular resonance lines are often blended with saturated ISM components, but lines from the 
excited levels do not suffer from this effect. The column densities of the excited levels can 
thus be used to determine the column density of the ground state, under the usual assump- 
tion of statistical equilibrium. Since the nebular components of the resonance transitions of 
S III at AA1012.50, 1190.20 are lost in the saturated ISM components and blends with other 
species, we use absorption lines of S III* and S III** to extrapolate the ground state column 
density, and subsequently derive the total ionic column density. The populations of these 
levels are also sensitive to n e , and we use the column density ratios to derive the line of sight 
electron density. 

Voigt component fits to S II and S III profiles are depicted in Figure 2a, with the fit 
parameters detailed in Table 1. We detected absorption from both excited levels of S III, at 
A1015.78 and AA1021.11, 1021.32, in the FUSE spectrum. Since only one S III* line provided 
useful information (the other lines in the FUSE and STIS wavebands are too severely blended 
and/or saturated to significantly constrain the fit), we chose to simultaneously fit S III* and 
S III**, fixing < v > and b, but allowing the column densities of the two levels to differ. 

The ground term of S II consists of a single level, so its ionic column density is given 
by the simultaneous fit of S II AA1250.58, 1253.81. S II A1259.52 is more heavily saturated, 
and does not significantly constrain the ionic column density when simultaneously fit with 
the other lines of the multiplet. The two lines we fit are mildly saturated, and therefore 
the derived S II column density is not well determined, but this ion does not contribute a 
significant amount of the total column of S. 



4-1.1. Electron Density from S III Level Populations 

In order to extrapolate the ground state population of S III from those of the excited 
levels, we used the 5-level task "nebular. ionic" by Shaw & Dufour (1995) in IRAF 3 , as- 
suming the electron temperature T e = 10, 500 K and density n e = 31, 600 cm -3 derived 
from nebular emission lines (DM01). However, at such a high density, the predicted ratio 
R 12 = N(S m*)/N(S III**) = 0.96 (hereafter, R {] refers to the population ratio of the i th 
to the j th levels above the ground state) is nearly a factor of two smaller than the observed 
value of 1.74±0.38. Note that saturation effects in S III* A1015.78 cannot explain the dis- 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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crepancy, for if this line were saturated, then i?i 2 would be even larger than the observed 
value. Furthermore, line-fit uncertainties are unable to bring the observed and predicted line 
ratios into accord - forcing the column densities to have a ratio of 0.96 leads to an inferior 
fit (higher x 2 ) for any value of b. We therefore explored the sensitivity of R12 to the adopted 
atomic data and physical parameters. 

Using more recent transition probabilities for S III (Nahar 1993) than the default values 
for the nebular. ionic task (Mendoza & Zeippen 1982) has almost no effect on the relative 
populations of the ground term levels. However, it should be noted that both sets of values 
assume LS coupling, which predicts the same relative oscillator strengths of lines within a 
multiplet. We also found that varying the electron temperature within a reasonable range 
(5-15 xlO 3 K) has little effect on the predicted S III level populations. On the other hand, 
R12 is sensitive to the assumed electron density, and increases as n c decreases. In Table 2, we 
show the predicted values of R m and R\ 2 for T e = 10, 500 K and a range of electron densities. 
The model agrees best with the observed R12 = 1.74 ± 0.38 at n c = 8800 cm~ 3 , with the 
l-o" errors allowing for the range n e = 6400-13600 cm~ 3 . It is clear that the observed S III 
column density ratio Ru is inconsistent with an electron density as high as reported by 
DM01, n c = 31,600 cm" 3 . 

It is important to note that DM01 derived the electron density from the optical doublet 
ratio [S II] A4072/A6725. The A4072 lines have high critical densities, and are therefore biased 
toward dense regions of the nebula. Other density diagnostics, such as [S II] A6717/A6731 
and [O II] A3726/A3729, imply n c > 10, 000 cm -3 , but are unreliable since they are near their 
high density limits. Therefore, the [S II] A4072/A6725 density indicator is the most reliable 
one measured by DM01. Using the same diagnostic ratio, dFP87 also derived a high electron 
density, n c = 10 5 cm~ 3 . It is possible that the more highly ionized portion of the nebula 
has a lower density than that occupied by singly ionized species. However, this argument is 
weakened by the fact that Flower et al. (1984) derived n c = 10 5 cm -3 from IUE observations 
of C III] AA1907, 1909. 

Is it possible that the nebular gas along the line of sight to the central star of SwSt 1 has 
a significantly lower density than other regions of the nebula? Since the absorption occurs 
along a "pencil beam" sample of the nebular gas in front of the central star, and thus samples 
a much smaller portion of the nebula than emission line studies, the electron density could 
be this low if the nebula is inhomogeneous. Other results (§§4.2, 5.1) indicate that this is 
indeed the case. 

We adopt n c = 8800±Hgg cm" 3 and T e = 10,500 K to derive #02 = 0.77 (see Table 2). 
This yields JV(S III A1012.50) = 0.77iV(S III**) = (6.41 ±2.32) x 10 14 cm" 2 , where the errors 
include the uncertainties in n c . The adopted column densities for the ground term levels of 
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S III are summarized in Table 3. 



4.2. Gas Phase Abundance of Iron 

Since the FUSE observations do not provide a means of directly determining the overall 
column density of ionized gas iV(H + ), we derive the Fe abundance in SwSt 1 by using sulfur 
as a reference element. S is particularly useful for this purpose, since it is not significantly 
depleted into dust in the diffuse ISM (e.g. Sembach & Savage 1996; Welty et al. 1999), nor 
is it affected by nucleosynthesis in the progenitor star. 

Because we do not observe all ionization stages of Fe, we use appropriate ionic ratios as 
surrogates for the elemental abundances. Fe II and Fe III have a similar ionization potential 
range to S II and S III (7.9-30.7 eV for Fe+ and Fe++, 10.4-34.8 eV for S+ and S++), and 
therefore we assume that iV(Fe + + Fe ++ )/iV(S + + S ++ ) ~ (Fe/S). We use the recommended 
solar abundances of Asplund et al. (2004) to derive the logarithmic abundance ratio relative 
to solar, [Fe/S] = log(Fe/S)-log(Fe/S) . 

As is the case for S III, the ground terms of Fe II and Fe III have multiple fine structure 
levels that may be significantly populated under nebular conditions. The column densities 
of some of these levels could not be measured directly, because absorption lines from the 
levels were undetected or saturated. Therefore, we computed the populations for these 
unobserved levels with models. We describe the methods used in the following subsections, 
before presenting the derived Fe abundance. 

4-2.1. Fe II Level Populations 

Only resonance transitions of Fe II are detected, despite the fact that this ion has several 
low-energy levels which may be significantly populated. The populations of the excited levels 
must therefore be predicted by models in order to correctly compute the total ionic column 
density. 

Iron is not included in the nebular. ionic task, so we use XSTAR (Kallman & Bautista 
2001) photoionization models to calculate Fe ionic level populations. XSTAR utilizes atomic 
data from the Iron Project (Hummer et al. 1993; see Bautista & Kallman 2001 and references 
therein), and computes populations of numerous levels for each ion. We assume a uniform 
temperature and density nebula, with a blackbody central star of T cff = 40, 000 K, and adopt 
other model input parameters from Table 12 of DM01. We tested the sensitivity of the level 
populations to the adopted parameters by varying elemental abundances, T e , radius of the 
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nebula, and central star temperature and luminosity. The ground term level populations of 
iron and sulfur ions are negligibly affected by these parameters, as is the derived abundance 
[Fe/S]. 

We computed models for several values of n c . From the S III fine structure levels 
(§4.1.1), we adopt n c = 8800 cm -3 , with an allowed range n e = 6400-13600 cm~ 3 . While 
Fe + may reside in the photodissociation region as well as the H II region, we find a relatively 
small total column density of Fe + compared to Fe ++ , so that uncertainties in the Fe II level 
populations have little effect on the derived Fe abundance. The resulting Fe level population 
ratios are listed in Table 3, along with the predicted column densities for unobserved or 
saturated lines. 

We find that a significant fraction of Fe + resides in excited levels for all of the considered 
n e (4000-10 5 cm" 3 ). The predicted column densities of the excited levels of the Fe II ground 
term are consistent with the measured upper limits of lines in the FUSE bandpass (Table 1). 
To compute the total column density of Fe + , we sum the calculated populations of the 
thirteen levels of the 6 D, 4 F, and 4 -D terms (higher energy levels contribute less than 1% to 
the Fe + abundance). The adopted column density for each level is given in Table 3. 



4-2.2. Fe III Level Populations 

Excited levels in the b D ground term of Fe III can give rise to several detectable absorp- 
tion lines in the FUSE bandpass. Fe III has no excited terms with energies below 2.4 eV, so 
it is sufficient to consider only the 5 D levels when computing the ionic column density. We 
detect lines arising from all ground term levels (Table 1), although the sole line of Fe in**** 
at Al 130.40 is quite weak. The column densities derived from the lines of Fe III and Fe HI* 
are very uncertain due to saturation, particularly that of the resonance line A1122.52. Since 
a large fraction of Fe III ions may reside in the two lowest energy levels, we use XSTAR 
and the absorption lines from the three highest energy 5 D levels to estimate their column 
densities. 

The predicted level population ratios are given in Table 3; in general, they correspond 
poorly to the observed ones. The discrepancies in Rq 2 and R12 are due to saturation of the 
resonance and Fe III* lines - the fitted column densities are lower limits. Therefore, we adopt 
the model column densities for Fe III and Fe III*. While -R32 is in good agreement with the 
model predictions, i? 42 is too small by a factor of 3-5. Although Fe in**** Al 130.40 is only 
marginally detected, it is very unlikely to be uncertain by such a large amount. Since FL42 is 
not sensitive to the adopted physical parameters of the photoionization models, we examine 
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the atomic data as a source of uncertainty. 

It is possible that the oscillator strengths, or /-values, of the Fe III a b D — z 5 P° multiplet 
are in error. While the Fe III /-values of Nahar & Pradhan (1996, hereafter NP96) are 44% 
larger than those recommended by Morton (2003; Ekberg 1993), the relative /-values within 
the multiplet are the same in the two studies, since both assume LS coupling and ignore 
relativistic corrections. If either of these assumptions does not hold, the relative /-values of 
lines within the multiplet would change. Alternatively, the adopted transition probabilities 
(NP96) and collision strengths (Zhang 1996) used by XSTAR to derive the level populations 
could be uncertain. However, the ground term level population ratios derived with this data 
(e.g. Keenan et al. 2001; Keenan 2004, private communication) agree to within 20% of those 
calculated with older atomic data (Keenan et al. 1992), showing that our results are not 
strongly dependent upon the adopted atomic data set. 

We are unable to explain the large discrepancy between the observed and predicted 
although Fe III Al 130.40 may have an uncertain /-value. With this caveat in mind, we use 
the observed i? 32 and R42 (which are from unsaturated lines), and the calculated values for 
R 02 and i? i2 . The adopted Fe III column densities are listed in Table 3. 

4-2.3. The Gas Phase Iron Abundance 

We compute the gaseous iron abundance using the sulfur and iron column densities 
summarized in Table 3, assuming that Fe/S « N(Fe + + Fe ++ )/N(S + + S ++ ). The gaseous 
Fe abundance is quite high, [Fe/S] = — 0.35±0.12, which implies that only about half of the 
nebular iron is depleted into dust along the line of sight to the central star. If one uses the 
Fe III /-values calculated by NP96 rather than those compiled by Morton (2003), [Fe/S] 
decreases to — 0.50±0.12. 

The derived [Fe/S] is shown in Table 4 for the adopted n e = 8800 cm -3 and the 1- 
o~ limits, 6400 and 13600 cm -3 ; note that the uncertainties in Table 4 are not the formal 
errors, since they do not include the uncertainty in n e . For comparison, we also show the 
Fe abundance derived with the n e value of DM01 (n c = 31,600 cm -3 ), and of dFP87 and 
Flower et al. (1984; n c = 10 5 cm -3 ). The inferred value of [Fe/S] decreases as n c increases, 
primarily because the Fe III ground state population relative to excited levels declines as gas 
density rises. Furthermore, if we have underestimated the saturation of the S II lines, the 
true value of [Fe/S] would be lower (because the actual N(S + ) would be larger than adopted 
here). Even given these considerations, the gaseous iron abundance along the line of sight 
is clearly high compared to the [Fe/H] observed in disk (—2.2 to —1.2 dex) and even halo 
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(—0.65 dex) diffuse interstellar clouds (Sembach & Savage 1996; Welty et al. 1999). 

We stress that the modest level of iron depletion derived from the FUSE spectrum is 
not necessarily characteristic of the whole nebula. In the following subsection, we show that 
iron is more heavily depleted in other regions of SwSt 1, by deriving the Fe abundance from 
optical emission lines. 

4-2.4- Comparison to Optical Emission Line Results 

From five optical [Fe III] lines, dFP87 computed [Fe/H] = -1.56 and [Fe/S] = -0.88, 
using their value of the S abundance. However, this determination relies on dated atomic 
data, and therefore may not be reliable. We re-derive the Fe abundance from optical emission 
lines, using current atomic data and the optical echelle spectrum of DM01, which was ob- 
tained with the Anglo-Australian Telescope (courtesy of O. De Marco; see also Appendix B 
of DM01). We measured line fluxes in the optical spectrum with the "splot" routine in 
IRAF, and de-reddened them using E(B — V)= 0.46 (DM01). The resulting line intensities 
relative to H[3 are reported in Table 5. 

We identify several [Fe II] and [Fe III] lines, but none from higher ions. Since most of 
the optical [Fe II] emission lines may be severely affected by fluorescence (Rodriguez 1999), 
we derive the Fe + abundance using only [Fe II] A8617.0, which is insensitive to fluorescent 
pumping (Lucy 1995; Baldwin et al. 1996). We identify 12 [Fe III] emission lines, which are 
of varying utility (see below) in determining the Fe ++ abundance. 

Several [Fe III] emission line ratios may be used as nebular diagnostics, as previously 
discussed by Garstang et al. (1978) and Keenan et al. (1993). However, the density di- 
agnostics are often degenerate, allowing for more than one value of n c (Bautista & Prad- 
han 1998; Zhang & Liu 2002). Using the diagnostic diagrams displayed in Figure 7 of 
Bautista & Pradhan (1998), we find that the ratios J(A4769.4)//(A4733.9) = 0.7 ± 0.1 and 
J(A4733.9)//(A4701.5) = 0.5 ± 0.1 indicate either logK/cm" 3 ) = 4.6 or 5.7. On the other 
hand, 7(A4754.7)//(A4733.9) = 0.9 ± 0.1 agrees with the other two ratios only for the low 
log(n c /cm~ 3 ) value; the high value significantly exceeds 6.0. Because of the agreement of all 
the [Fe III] density diagnostics at log(n c /cm~ 3 ) = 4.6, and its close concordance to the value 
derived by DM01, we adopt this value for the electron density. Unfortunately, none of the 
temperature-sensitive transitions discussed by Zhang & Liu (2002) were detected. Therefore 
we adopt T e = 10, 500 ± 500 K (DM01) to derive the Fe ionic abundances. 

In order to calculate the Fe + and Fe ++ abundances (Table 5), we use XSTAR models 
(see §4.2.1) to derive the fractional population in the upper level of each transition. Note 
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that some of the [Fe III] lines give discrepant ionic abundances compared to other transitions 
(Table 5). In particular, [Fe III] AA4607.0, 4777.7, 4881.0, 4930.5, and 4987.9 indicate signif- 
icantly different ionic abundances than the other seven lines. Rodriguez (2002) noted that 
AA4607.0 and 4777.7, which yield higher Fe ++ abundances, may be contaminated by blends 
with other species; in addition, [Fe III] A4930.5 may be contaminated by [O III] A4931.2. 
However, blends cannot explain the anomalously low fluxes of A4881.0 and A4987.9. Mea- 
surement uncertainties are unlikely to be the cause for the inconsistency, particularly in the 
case of the relatively strong line A4881.0. Interestingly, both of these lines originate in the 
a 3 if 4 level. 

In Table 5, we compare the observed line intensities to those predicted by Keenan et 
al. (2001; K01) and Keenan et al. (1992; K92, pre-Iron Project) to examine the effects of 
different sets of atomic data. While the predicted line intensities are similar in most cases, 
the K92 values actually reproduce the a 3 H — 5 D intensities quite well. By comparison, K01 
(who use the same atomic data as we do here) predict /(A4881.0) and /(A4987.9) to be 
too large by a factor of two or more. It is possible that the atomic data for the a 3 H — 5 D 
multiplet is not reliable, but it should be noted that neither Rodriguez (2002) nor Perinotto 
et al. (1999) commented on inconsistencies from these two lines. In any case, we do not use 
AA4607.0, 4777.7, 4881.0, 4930.5, and 4987.9 in the derivation of the Fe++ abundance. 

We derive Fe ionic abundances (Table 6) using [Fe II] A8617.0 and the seven most reliable 
[Fe III] lines, assuming log(n c /cm~ 3 ) = 4.6 and T e = 10,500 K. We vary log(n e /cm~ 3 ) from 
4.5 to 5.0, in accord with previously derived electron densities (DM01; dFP87; Flower et 
al. 1984), and assume an uncertainty of ±500 K for T e , as given by DM01, to produce the 
error bars on the ionic abundances listed in Table 6. As discussed above, Fe + is a trace ion 
in the nebula, so the derived Fe abundance is not strongly affected even if much of the Fe + 
resides in the photodissociation region, under significantly different physical conditions than 
we assume. 

In order to correct for unseen stages of ionization (primarily Fe +3 ), we adopt the ion- 
ization correction factor (ICF) given by Rodriguez (2002): 

Fe/H = 1.1 x (0/0+) x (Fe ++ /H + + Fe + /H+). 

Using data from Tables 10 and 11 of DM01 (we adopt the "empirical abundance" for O in 
their Table 11), we find ICF(Fe) = 1.1 x (0/0 + ) = 1.16 ±0.56. The relatively large error is 
due to the uncertainty in the total O abundance. This results in [Fe/H]= —1.64 ± 0.24, or 
[Fe/S]= —1.15 ±0.33, using the S abundance from DM01. In other words, the Fe abundance 
we derive from the optical data of DM01 agrees well with that of dFP87. Iron is depleted by 
a factor of 10 or more according to the optical data, which starkly contrasts with the very 
modest depletion we find from the FUSE data. 
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Thus, it appears that the high gaseous Fe abundance that characterizes the column 
toward the central star is not representative of the emitting ionized gas. This result is 
strengthened by the fact that the [Fe/S] derived from the FUSE spectrum is not strongly 
dependent upon the assumed n c (Table 4). 

4.3. Abundance of the s-process Element Ge 

The slight nebular Fe depletion along the line of sight implies that other refractory 
elements will also be relatively undepleted. SDB derived the Ge abundance in four PNe 
(including SwSt 1) from the resonance line Ge III A1088.46 in the FUSE band, and reported 
a range of possible Ge abundances, depending on the unknown degree of depletion. The 
range found for SwSt 1 spanned [Ge/S] = (0.86-1.48)±0.15 (note the typographical error 
in Table 2 of SDB), depending on whether Ge is wholly gaseous, or is depleted to the same 
extent as in cold interstellar clouds (—0.62 dex; Savage & Sembach 1996). The modest Fe 
depletion indicates that [Ge/S] is at the low end of the range found by SDB. 

We have re-evaluated the Ge abundance with our new method of analysis (§3.1) and 
updated atomic data from Morton (2003). Using the ionization correction factor Ge/S = 
iV(Ge ++ )/iV(S + + S ++ ), based on similarities of the ionization potential ranges, we find 
[Ge/S] = 0.72±0.06. This agrees within the errors with the value [Ge/S] = 0.86±0.15 
calculated by SDB. 

Many theoretical and observational studies of the s-process discuss n-capture element 
yields relative to iron. Assuming that Ge/Fe = iV(Ge + +)/iV(Fe ++ ), we derive [Ge/Fe] = 1.08± 
0.13 using Fe III /-values from Morton (2003), or 1.24±0.13 using values from NP96. (Note 
that [Ge/Fe] is not equivalent to [Ge/S] - [Fe/S]; the ionization potential range of Fe + does 
not coincide with that of Ge ++ , and is therefore not used to derive [Ge/Fe]). In Table 4, we 
display the Ge abundance relative to S and Fe for the range of densities implied by UV S III 
absorption lines, as well as the n e found from optical emission lines. If iron is depleted, then 
the derived value of [Ge/Fe] is an upper limit. However, [Ge/S] is enhanced by at least a 
factor of 5 for the entire range of electron densities we consider. 
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5. RESULTS FOR THE NEUTRAL/MOLECULAR REGION 

5.1. Non-Detection of Nebular H 2 

The original motivation for FUSE GO program B069, for which this observation of 
SwSt 1 was taken, was to study the molecular components of PNe through examination of 
the H 2 Lyman and Werner resonance lines seen in absorption against the continua of their 
central stars. Accordingly, targets were selected to have velocities well-separated from low- 
velocity foreground interstellar material, and to display near-infrared quadrupole emission 
lines from vibrationally-excited H 2 or some other indicator of neutral or molecular material. 
SwSt 1 was selected as a target for the H 2 search because it meets both of these criteria. 

5.1.1. Previous Evidence for Neutral/ Molecular Material in SwSt 1 

SwSt 1 was found by Dinerstein et al. (1995b; hereafter DSU) to possess nebular com- 
ponents in the Na I "D" lines at 5890.0, 5895.9 A. Since neutral Na has a very low ionization 
potential (5.1 eV) and appears to correlate with other indicators of atomic and molecular 
material (see Table 5 of DSU), Na I is presumed to be a tracer of neutral material in nebulae. 

In the case of pure absorption by foreground material, the Na I lines can in principle be 
used to infer a column density of circumstellar material from the classical "doublet ratio" 
method (e.g. Dinerstein & Sneden 1988), which has been widely applied to diffuse interstellar 
clouds. However, this simple approach is probably invalid for SwSt 1, where the modest 
equivalent widths of the Na I lines (55±15 mA) appear inconsistent with their ratio, which 
is close to the fully saturated value of unity. The reason for this inconsistency is evident from 
Figure 5 of DSU, where it can be seen that the nebular absorption in SwSt 1 is partially 
filled in by scattered Na I photons; the true absorption profiles are actually deeper than they 
appear. The observed absorption minima fall at a heliocentric radial velocity of -47 km s -1 , 
blueshifted from the systemic emission line value of -20 km s _1 ; the blending of the emission 
and absorption components is not likely to affect the apparent velocity by more than a few 
km s _1 . 

Circumnebular H I 21 cm absorption was subsequently reported by Gussie & Taylor 
(1995), at a velocity of ~-38 km s" 1 . They did not, however, detect CO J = 2-1 millimeter 
emission, although all of the other four PNe they detected in H I were also seen in CO. 
Huggins et al. (1996) also searched for CO J = 2-1 in SwSt 1 and did not detect it. This 
suggests that SwSt 1 possesses only a modest amount of neutral/molecular material. 

SwSt 1 is a source of near-infrared H 2 emission, although its spectrum has not been 
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studied as intensively as some other PNe. In their Appendix B, DM01 reported a detection 
of the 2.122 /iin v = 1-0 S(l) line at a strength about 5% that of H I Br7. Additional H 2 lines, 
in particular v = 1-0 S(0) 2.223 jum and several Q-branch lines, are visible in the spectrum 
of DM01, although the line fluxes are not cited in their Table B2. Likkel et al. (2004) 
also detect 2.122 and 2.223 /xm, but no other H2 lines; their data do not extend longward 
of 2.26 /iin and therefore do not include the Q-branch. They measure a 2.223/2.122 /im 
flux ratio consistent with radiative excitation of moderately dense gas in a strong radiation 
field. If radiative excitation is occurring, the Lyman and Werner lines in the FUSE band, 
which initiate the fluorescent cascades, should be detectable as long as there is H 2 along the 
absorption column. 



5.1.2. Upper Limits on N(H 2 ) for Specific J Levels 

We do not detect nebular H 2 absorption in the FUSE spectrum of SwSt 1. In order to 
quantify this statement, we measure limits on the equivalent widths of nebular components 
of specific H 2 transitions on the blue wings of the corresponding interstellar features. Because 
the velocity shift is modest (only 2-3 spectral resolution elements), this is not feasible for lines 
that are very saturated in the ISM, namely those arising from the J = and 1 levels. Instead, 
we examine lines from J > 2, which have less saturated and hence narrower interstellar 
components. We select several lines, arising from each of the J = 2, 3, 4, and 5 levels, that 
are free from interference by other species and fall in regions with good S/N ratio in the 
continuum. 

The selected transitions are listed in Table 7, and we display a representative sample of 
these lines in Figure 3. For comparison, we show J = and 1 lines in the upper panels, from 
which it can be seen that blending with the interstellar component (thick solid line) prevents 
meaningful measurements of limits on the nebular column density. We use two methods 
to set upper limits on the equivalent widths of the nebular components. In some cases, we 
were able to "fit" a nebular feature by interpolating a conservatively high continuum near 
the nebular velocity. For an alternate estimate, we evaluated the RMS noise level in the 
adjacent continuum, assumed a line width equal to the spectral resolution element, and took 
the corresponding 3-<r value as our upper limit. The W\ values from both methods, typically 
~ 5-15 mA, were converted to column densities iVj by assuming that the lines are on the 
linear portion of the curve of growth. The Nj limits are given in Columns 4 and 5 of Table 
7, and depicted in Figure 3 by dotted lines showing Voigt profiles with the measured iVj 
limit at the nebular velocity, assuming b = 11.4 km s" 1 (the value measured for nebular O I 
lines). Overall, our upper limits are approximately JVj ~5x 10 13 cm" 2 for J = 2 and 3, and 
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JVj ~4 x 10 13 cm" 2 for J = 4 and 5. 

5.1.3. Upper Limits on the Total i/ 2 Column Density 

In order to translate these level-specific column densities into a total H 2 column, it is 
necessary to make assumptions about the excitation equilibrium. Given fj, the fraction of H 2 
molecules in level J, then it is possible to compute N tot = Nj/fj for each J. For purposes 
of a rough estimate, we assume that the lowest few J levels are dominated by collisional 
processes, so that their populations are predicted by Boltzmann statistics. To normalize 
the level population Nj = gj e - {xj/kT \ we need the appropriate partition function Z(T) = 
£ (iVj) over all J levels that contain significant populations. 

We compute partition functions for H 2 for several representative values of the rotation 
temperature T rot : 200, 500, 800, and 1600 K. We include the 20 energy levels with the largest 
populations, which consist of the J = - 15 levels of the v = vibrational state, and J = 
0-5 for v — 1, which begins to interleave the energy level ladder of v — between J = 8 
and 9. For radiative excitation, a single value of T rot does not strictly apply; however, the 
effective values of T rot for J < 5 for moderately high radiation field and gas densities lie in the 
range ~ 200 - 500 K (Draine 2004, private communication). Furthermore, in BD+30°3639, 
which has highly excited circumstellar H 2 , T rot ~ 800 K for v = 0, while the v > 2 levels 
have substantially higher values of T rot (Dinerstein & Bowers 2004; Dinerstein et al. 2005, in 
preparation). 

In Table 8 we convert the upper limits on the specific J levels into the corresponding 
total column densities iV tot (H 2 ). In general, the strongest constraints on the total amount 
of H 2 come from the J levels with the largest fractional populations, but the J value with 
the largest fraction varies with temperature. For T rot > 500 K the tightest constraint comes 
from the J = 3 lines: A^ tot (H 2 ) ~2 x 10 14 cm -2 . For the lowest temperature we calculate 
here (more appropriate for lower-density gas and a weak radiation field), the best limit, 
N to t(H 2 ) < 7 x 10 14 cm" 2 , is from J = 2. 

It should be noted that the derived upper limit to A^ t(H 2 ) may be in error if the H 2 
is sufficiently cold that only the J = and 1 states are significantly populated, since we 
were not able to derive meaningful upper limits to the column densities of these levels. For 
example, if T rot = 100 K, then 1.2% of the H 2 would be in the J = 2 state, so that our upper 
limit of iVtot(H 2 ) = A^ 2 // 2 < 7 x 10 14 cm" 2 (derived for T rot = 200 K) would be too low by 
a factor of about 6. Temperatures lower than 60 K would be necessary for our estimate to 
be too low by two orders of magnitude. However, given that H 2 in the diffuse ISM, far from 
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strong UV radiation fields, exhibits rotational temperatures of >80-90 K (Shull et al. 2004; 
Savage et al. 1977), it is very unlikely that circumstellar material close to the luminous, hot 
central star of a PN would be this cold. 



5.1.4- Comparison to H 2 Seen in Emission 

The flux of H 2 1-0 S(l) 2.122 /im was reported as (1.5 ± 0.2) x 10~ 13 erg cm" 2 s" 1 by 
DM01, and (1.1 ±0.3) x 10" 13 erg cm" 2 s" 1 by Likkel et al. (2004). While there are no direct 
determinations of the size of the H2-emitting region, the similarity of the fluxes measured in 
different aperture sizes implies that it must be relatively compact. HST images of the ionized 
gas suggest a knotty ring of approximate dimensions 1'.'3 x 0'.'9, or 4.3 x 10~ n ster (DM01). 
Since in most PNe the H2-emitting regions are slightly more extended than the ionized gas 
(e.g. Latter et al. 1995), we assume that the H 2 fills solid angle ~5 x 10~ n ster. This corre- 
sponds to a surface brightness for the 2.122 /im line of (2.6 ± 0.4) x 10~ 3 erg cm~ 2 s _1 ster -1 . 
Even this large value, 20 times the surface brightness of the "canonical" fluorescent H 2 PN 
Hubble 12 (Dinerstein et al. 1988; Ramsay et al. 1993), is only a lower limit, since it is 
beam-averaged and the emission may not be uniform. 

For optically thin emission, this surface brightness implies a beam-averaged column 
density in the originating level of the transition (y — 1, J = 3) of 1.0 x 10 17 cm" 2 . Even in 
a highly excited cloud, this level contains only a small fraction of the total H 2 : for thermal 
excitation, f 1:3 (800K) ~ 2xl0~ 4 and /i, 3 (1600iir) ~ 7xl0" 3 . Therefore the total column 
density of H 2 would have to be at least 7V tot (H 2 ) > 10 19 cm -2 . Comparing this estimate 
to our limits on the line of sight iV tot (H 2 ) from the UV data, we see that they differ by at 
least four orders of magnitude. If the absorption line of sight contained the column of H 2 
indicated by the infrared emission lines, the nebular components of the UV lines would be 
deeply saturated, even if the H 2 is very hot (T rot = 1600 K; dashed lines in Figure 3). This 
disparity indicates an extreme asymmetry or dumpiness in the distribution of molecular 
material in SwSt 1. 



5.2. Excitation of O I in SwSt 1 

5.2.1. UV Absorption Lines from Excited O I 

Despite the absence of nebular H 2 absorption in SwSt 1, some neutral species are clearly 
present in the FUSE spectrum. Some of the low-ionization species we detect, such as Ar I, 
Fe II, and S II, may exist in both the ionized nebula and in the transition region between 
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ionized and neutral material (the photodissociation region, or PDR). However, we also detect 
absorption from O I, which has an ionization potential nearly identical to that of H I, and 
therefore is expected to reside primarily in the neutral zone. 

Under the conditions of a PDR, the excited fine structure levels of the O I 3 P ground 
term may be significantly populated, and therefore can be detected in absorption along with 
resonance transitions. If the excitation of these levels is exclusively due to collisions, one 
can use their column density ratios as a measure of the kinetic gas temperature. While the 
lines of the well-known resonance triplet at 1302-1306 A in the STIS waveband suffer from 
saturation and blending, the FUSE band offers a variety of alternate O I multiplets, with 
a range of different /-values and hence sensitivities to saturation effects. Many of these lie 
below 1000 A, where the spectrum is crowded and has low S/N, so we focus on the multiplet 
near 1040 A. Although the ground state line (A1039.23) cannot be used due to saturation 
and blending with the ISM component, simultaneous fits show that the absorption lines from 
the excited levels in this multiplet are free from such effects (see Figure 2a). 

Taking the ratio of the column densities derived from these fits (Table 1), we find 
iV(0 T* A1041.69)/7V(O I* A1040.94) = 0.71±0.16. This value exceeds by about a factor of 
two the ratio of statistical weights (1/3), which is the limiting value for LTE in the infinite 
temperature case, T >> Ax/k (where Ax is the energy difference between the two levels). 
This population inversion demonstrates incontrovertibly that collisions are not the primary 
excitation mechanism for the fine structure levels of O I in SwSt 1. 



5.2.2. Alternate Excitation Mechanisms 

Several processes other than collisional excitation can populate the excited levels of O I. 
These processes were discussed and evaluated in detail by Grandi (1975a, b), in order to 
understand the strength of O I A8446.5 and other permitted lines of O I in the Orion nebula. 
These excitation mechanisms include recombination into excited levels; fluorescent excitation 
by the stellar continuum; or fluorescent excitation by H I Lyman lines, primarily Ly (3 but 
also to some extent Ly e and Ly (. In order to distinguish among these possibilities, we 
examine the intensities of various optical and near-infrared O I lines in SwSt 1 as reported 
in the literature. 

We use the data of DM01 (see §4.2.4 for more details) to measure optical and near- 
infrared O I lines. The resulting ratios of various O I line fluxes to A8446.5 are summarized 
in Table 9. We also list the corresponding ratios in Orion (Esteban et al. 2004), for which 
Grandi (1975a, b) concluded that fluorescent excitation by the stellar continuum was the 
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only mechanism that could explain the observed pattern of lines. 

Recombination can be ruled out as the dominant O I excitation mechanism in SwSt 1 
by virtue of the strength of O I A8446.5 relative to Ha: the line flux ratio is ~0.003, over a 
factor of 30 larger than predicted for recombination (Rudy et al. 1989). Additional support 
comes from the weakness of O I A7774.2 relative to A8446.5; if recombination populated these 
levels, these two lines should have comparable strengths. Starlight excitation is usually cited 
as the mechanism which produces O I A8446.5 in Galactic nebulae (Grandi 1976), but Rudy 
et al. (1989) have argued, on the basis of measurements of near-infrared O I lines, that Ly (3 
fluorescence is responsible for exciting O I in the compact PN IC 4997. However, Lyman 
line fluorescence cannot explain the strengths of O I lines from states other than 3 D, such as 
AA5554.9, 6046.4, and 13164 (Grandi 1975a; Rudy et al. 1989), which cannot be populated 
via line excitation. 

The general consistency between the O I line intensities in SwSt 1 and Orion seen in 
Table 9 strongly supports the conclusion that the excitation of the higher levels of O I is 
dominated by stellar continuum fluorescence. Since the atoms will cascade to lower levels, 
eventually reaching the fine structure levels of the ground 3 P term, we suggest that the 
population inversion we find in the FUSE spectrum is a result of this fluorescent cascade 
initiated by stellar continuum photons. 

5.2.3. Implications for the Interpretation of the Far-Infrared [O I] Emission Lines in PNe 

[O I] 3 Pi - 3 P 2 63 /iin, the transition from the first excited fine structure level to the 
ground state, is one of the strongest cooling lines in PDRs. Although [C II] 158 fim is the 
strongest line at low densities, the 63 /iin line dominates the cooling of PDRs with intense 
incident UV fields and high densities, conditions that apply to the circumstellar PDRs of 
PNe (Tielens & Hollenbach 1985; Hollenbach & Tielens 1997). The companion [O I] 3 Po ~ 
3 Pi line at 145 /im is weaker than 63 fim due to the higher energy and lower statistical 
weight of its originating level. Under the assumption that collisional excitation dominates 
the level populations, ratios among the three lines can be used as diagnostics of the local 
thermal parameters n and T (e.g. Watson 1984). 

However, the populations in the O I ground term can be redistributed by fluorescent 
excitation in a strong stellar radiation field. As seen above, stellar continuum fluorescence 
can overpopulate the upper level of the [O I] 145 /im line relative to its value from purely 
collisional excitation. Therefore, when starlight excitation affects the ground term levels of 
O I, the true PDR temperature will be lower than indicated by the line flux ratios of [O I] 63, 
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145, and [C II] 158 fim. In the gas column of SwSt 1 probed by FUSE, starlight fluorescent 
excitation is so dominant that there is no ambiguity about its importance. However, in other 
PDRs, more modest contributions might be present but not easily recognized. 

Observations with NASA's Kuiper Airborne Observatory (KAO) and the Infrared Space 
Observatory (ISO) have been used to determine the physical parameters of PDRs for sev- 
eral PNe, although SwSt 1 itself has not been observed. In order to interpret the line 
ratios, a correction must be made for the elevated C/O abundances of these objects. Ap- 
plying such corrections to measurements from the KAO, Dinerstein et al. (1995a, 1997) 
found log(n/cm~ 3 ) ~ 4.3, T > 1000 K for the PDRs associated with BD+30°3639, IC 418, 
NGC 6572, and NGC 7027. (We cite only a lower limit for T, since the ratios become in- 
sensitive to the precise value for T > 1000 K.) It is interesting to note that O I optical and 
near-infrared line ratios indicate that stellar continuum fluorescence is important in all four 
of these PNe (BD+30°3639 - Aller & Hyung 1995, Rudy et al. 1991b; IC 418 - Sharpee 
et al. 2003; NGC 6572 - Hyung et al. 1994, Rudy et al. 1991a; NGC 7027 - Baluteau et 
al. 1995, Rudy et al. 1992). Liu et al. (2001) observed the same four objects as well as 
18 additional PNe, with the Long Wavelength Spectrometer on ISO. With the exception of 
NGC 7027, the measured flux ratios of Liu et al. differ by up to a factor of two from those 
found by Dinerstein et al., yielding lower temperatures (T ~ 200-400 K) and much higher 
densities (log(n/cm -3 ~ 5.3 — 5.5) for these PNe. The reason for this difference is not clear. 
One possible source of the discrepancy is the correction for C/O, which is uncertain by as 
much as a factor of 2-3 in PNe. Another factor is that the ISO observations were taken 
with a larger entrance aperture than the KAO data, and may include a greater contribution 
from more extended, cooler gas, where the stellar radiation field is weaker and collisional 
excitation dominates. Thus, the contribution of the fluorescently excited component may be 
diluted compared to that seen by KAO. 

This effect has probably been seen, albeit not recognized, in other sources such as 
the p Oph cloud (Liseau et al. 1999a) and Carina nebula (Mizutani et al. 2004), where 
excess emission is seen in the 145 /zm line compared to the predictions of PDR models and 
other diagnostics of gas temperature. Indeed, Liseau et al. (1999a) specifically suggest that 
an unidentified nonthermal pumping mechanism may be responsible (also see Liseau et al. 
1999b). Our FUSE observations of SwSt 1 provide an explanation for this widespread puzzle. 



-22 - 



6. DISCUSSION 

6.1. Nature of the Inhomogeneities in SwSt 1 

Many of the physical properties that we derive from the FUSE spectrum of SwSt 1 - the 
electron density, the high gaseous iron abundance, and the absence of nebular H 2 absorption 
- differ dramatically from values derived previously from emission lines. However, FUSE 
probes only a "pencil-beam" of gas in the line of sight to the central star, a small portion 
of the nebula compared to emission line studies; the differences therefore suggest that the 
absorption line of sight is not representative of conditions elsewhere in the nebula. These 
results indicate a non- uniform distribution of ionized and molecular material in SwSt 1, 
either due to a highly asymmetrical global structure, or extreme clumping. 

It is possible that global structure can explain the observations. The images of the 
optical emission in DM01 suggest a toroidal morphology, tipped with respect to our line of 
sight. (Assuming that the optical nebula is really a circular ring, the observed aspect ratio 
implies an inclination of ~60° from pole-on.) In this case, the FUSE spectrum looks through 
a polar cone in which the gas has mostly been evacuated. The line of sight to the central 
star, which traverses this cone, would then be a region deficient in both molecules and dust, 
where the gas density is low compared to the toroidal ring. 

Alternatively, small-scale clumping may be responsible for the nebular inhomogeneities. 
In fact, such dumpiness may be quite common in PNe. For example, O'Dell et al. (2002) find 
numerous knots in nearby PNe, and Gongalves et al. (2001) discuss 50 PNe with detected low 
ionization structures. Small-scale structure may form quite early in the evolution of a PN, 
from instabilities created during the interaction of the fast and slow winds (Dwarkadas & 
Balick 1998) or with the passage of the ionization front (e.g. Garcia-Segura & Franco 1996). 
Condensations in the progenitor's envelope during the AGB phase could also contribute to 
knots in the ejecta (Dyson et al. 1989), although Huggins & Mauron (2002) find a lack of 
evidence that this occurs. 

Recent studies have shown that PNe with H-deficient central stars (such as SwSt 1) may 
be particularly likely to form clumpy nebulae. Acker et al. (2002) and Peha et al. (2003) 
find that PNe with [WC] central stars tend to exhibit a higher degree of turbulence (which 
implies the presence of small-scale instabilities) than those with H-rich stars. Indeed, the high 
metal abundance and mass-loss rates of [WC] central stars likely prolong the momentum- 
conserving phase of the nebula (Acker et al. 2002), which is prone to nonlinear thin shell 
instabilities (Dwarkadas & Balick 1998). Furthermore, clumpy winds from [WC] central 
stars, as detected by Grosdidier et al. (2000), can create or enhance inhomogeneities in the 
nebula. While Grosdidier et al. did not detect inhomogeneous outflows from the central 
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star of SwSt 1, they acknowledge that their observations were insufficient to detect wind 
variability on timescales longer than about an hour. However, they detected clumpy winds 
from the central stars of NGC 40 and BD+30°3639, which are similar to that of SwSt I. 

Since the angular scale of SwSt 1 is so small that even HST images cannot fully resolve 
its intrinsic structure, it is not possible to assess whether the non-uniform distribution of 
nebular material is due to global asymmetries or clumping. However, it is clear that the 
nebula is inhomogeneous, and the sight line to the central star does not intersect material 
with the properties exhibited by emission lines. 

6.2. Inhomogeneities of the Dust and Molecular Gas in SwSt 1 

The spatial segregation of H 2 into dense regions, implied by the absence of UV H 2 
absorption in an object with near-IR H 2 emission, is not unexpected in view of the fragility 
of H 2 molecules to photodissociation and the intensity of the stellar UV radiation field. 
Tielens (1993) estimated that the lifetime of H 2 molecules outside of dense knots - where 
self shielding and dust extinction protect them from photodissociation - is on the order of a 
few years, and therefore predicted that H 2 in PNe must be confined to dense clumps. Recent 
high spatial resolution observations of H 2 and CO emission in the Helix and Ring Nebulae, 
where small linear dimensions can be resolved, have confirmed this expectation (Huggins et 
al. 2002; Speck et al. 2003). Even given the relative youth of SwSt 1 (estimated as ~ 300 
yr by DM01), the timescales for the propagation of photodissociation fronts are sufficiently 
short (e.g. Hollenbach & Natta 1995) that the H 2 will be largely dissociated along paths 
which do not contain condensations with N > 10 19 " 20 cm~ 2 . 

The origin of the inhomogeneities in the gas phase Fe abundance and the corresponding 
dust-to-gas ratio in SwSt 1 is less clear. Our results show a difference of ~ 1 dex between 
the gas phase Fe abundance in the UV absorption column and the emitting gas, where it is 
depleted by about a factor of ~ 40, similar to values found in the warm ionized ISM and H II 
regions (Sembach & Savage 1996; Welty et al. 1999; Rodriguez 2002). The low optical Fe ++ 
abundance is corroborated by Likkel et al. (2004), who find similarly low Fe ++ /H + from 
near-infrared [Fe III] lines in SwSt 1. Further evidence for the presence of dust in SwSt 1 
comes from the detection of thermal dust emission in the infrared. SwSt 1 exhibits the dual 
dust chemistry that is typical of PNe with [WC] central stars (Cohen et al. 2002), as indicated 
by the presence of features from both O-rich (Aitken et al. 1979; Zhang & Kwok 1990) and 
C-rich dust (Roche et al. 1996; Szczerba et al. 2001). The O-rich grains dominate by mass: 
silicates account for about 95% of the dust, with the other 5% giving rise to the aromatic 
emission bands usually attributed to PAHs (Casassus et al. 2001). On general considerations, 
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Fe should be strongly depleted in regions with abundant silicate grains (Draine 2003), and 
has also been found to be underabundant in nebulae with C-rich dust (e.g. Shields 1978; 
Perinotto et al. 1999). 

There are two possible scenarios that could produce the observed difference in Fe deple- 
tion between the absorption column and the emitting regions. If the ejected material initially 
had a uniform dust-to-gas ratio, then the dust must have been destroyed more efficiently in 
lower density regions of the nebula. Alternatively, variations in dust formation or growth 
efficiency in a clumpy outflow could have caused a non-uniform dust-to-gas ratio, and the 
initial inhomogeneity remains imprinted on the evolving nebula. 

The main dust destruction agents in nebulae are shocks and energetic stellar photons. 
Although the central star of SwSt 1 drives a fast wind (900 km s -1 ; DM01), the wind-shocked 
gas should lie in the central cavity, and be highly ionized and heated to T ~ 10 5-6 K (e.g. Chu 
et al. 2004). Shocks moving into the dense optically emitting material would be much slower 
(Tielens 1993), and therefore less able to destroy the dust. The fact that we measure radial 
velocities for the UV absorption lines that are consistent with the overall expansion velocity 
of the optical emitting gas (DM01) argues that the absorbing gas was not accelerated by a 
shock. In addition, the ionization balance seen from the UV absorption lines is consistent 
with that indicated by the emission lines. Therefore, based on kinematic and excitation 
arguments, shock processing of the gas and dust in SwSt 1 seems improbable. 

Unlike H 2 molecules, dust grains are relatively robust with respect to evaporation by 
stellar radiation, with the exception of small grains (containing fewer than 20-30 atoms) 
that have small heat capacities relative to the energies of typical stellar UV photons. These 
small grains cannot re-radiate the absorbed photon energy via an equilibrium process, and 
can undergo large temperature fluctuations. Such grains, whether carbon or silicate-based, 
can be heated to temperatures exceeding the sublimation temperature of refractory elements 
after absorbing just one or two high-energy photons (e.g. Guhathakurta & Draine 1989). 

Observational studies have drawn mixed conclusions about radiative destruction of dust 
in PNe. Observed gas phase depletions of elements such as Fe and Mg do not show any clear 
trends with nebular age (Middlemass 1988; Perinotto et al. 1999). Similarly, it does not 
appear that dust-to-gas ratios derived from infrared emission are lower in older, more evolved 
PNe (Hoare et al. 1992). Such arguments depend on comparing different objects, which may 
not be appropriate if the initial dust formation efficiency differs among their progenitor stars. 
Another approach is to search for evidence of dust destruction within specific nebulae. In this 
regard, the fact that the dust-to-gas ratios derived for the ionized and neutral (PDR) regions 
in NGC 7027 agree to within a factor of two is taken as evidence against dust destruction in 
the ionized gas by Hoare et al. (1992). On the other hand, there is evidence for a "depletion 
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gradient" in some PNe in the sense that some refractory elements (e.g. Al, Ca, and Mg) are 
less depleted in regions of higher than those of lower ionization, suggesting that atoms have 
been released from grains into the gas phase by the stellar radiation field (Kingdon et al. 
1995; Kingdon k Ferland 1997). 

In SwSt 1, we find very disparate gas phase Fe depletions in different regions of the 
nebula with the same ionization states, which makes it seem unlikely that there has been 
a substantial difference in the amount of radiative dust destruction between these two re- 
gions. This indicates that the disparate Fe depletions are most likely caused by initial 
inhomogeneities in the dust formation or growth efficiency in the AGB outflow. Such inho- 
mogeneities could be caused by low temperature structures (e.g. starspots) in the atmosphere 
of the progenitor AGB star, as suggested by Frank (1995) and Soker & Clayton (1999). In 
addition, the grain growth efficiency is higher in denser regions of a clumpy outflow (e.g. 
Dorfi & Hofner 1996), so that the grain size distribution would be skewed toward larger radii 
in dense gas. If the density inhomogeneities are preserved into the PN phase, dust in the 
more rarefied gas is more susceptible to photodestruction because of the smaller grain size 
distribution. This implies that radiative dust destruction could have been more efficient in 
lower density regions of SwSt 1, even though the incident radiation field is similar for the low 
and high density gas (as indicated by the comparable ionization balance of the absorbing 
and emitting gas). 

The dust deficit in the gas probed by the UV absorption lines may have a direct influence 
on the amount of H 2 present, since the primary channel for H 2 formation is through grain 
surface reactions. The removal of this mechanism for reforming H2 is likely to lower the equi- 
librium H 2 abundance that can be maintained in the presence of a dissociating and ionizing 
radiation field (Aleman & Gruenwald 2004 and references therein). Therefore, although we 
attribute the extreme contrast in the H 2 abundance primarily to photodestruction, it is pos- 
sible that the dust inhomogeneities amplify and maintain the spatial variations in molecular 
material. 



6.3. Nucleosynthesis and Dredge-Up in the Progenitor Star 

Germanium (Z = 32) isotopes can be produced by slow neutron (n) -capture nucleosyn- 
thesis (the s-process) in PN progenitor stars during the AGB phase of evolution. Nuclei 
created by the s-process may be convectively transported to the stellar envelope along with 
carbon during third dredge- up (TDU; e.g. Busso et al. 1999). The fact that we find an 
enhanced Ge abundance relative to solar in SwSt 1 implies that the s-process and TDU were 
efficient in the progenitor star. 
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We measure a Ge enrichment [Ge/S] = 0.72±0.06, or [Ge/Fe] = 1.08 ±0.13, which can 
be compared with theoretical models of stellar s-process enrichments, given the metallicity. 
The S abundance derived by DM01 and dFP87 implies that the progenitor star of SwSt 1 
may have had a metallicity as low as ~ 0.3Z & . However, the O and Ne abundances indicate 
that the metallicity is closer to solar. The derived Ge overabundance is slightly larger 
than predicted by recent theoretical models of the s-process in this metallicity range (Busso 
et al. 2001; Busso 2003, private communication; Goriely & Mowlavi 2000). However, many 
uncertainties exist in the strength of TDU and in the effect of convective overshoot and 
rotation on the s-process (e.g. Herwig et al. 2003; Siess et al. 2004). Therefore, we consider 
the derived Ge enrichment to be consistent with evolutionary models, within uncertainties 
in the metallicity and yields. 



7. CONCLUSIONS 

We now summarize the principal results of our study: 

1. We have analyzed the FUSE spectrum of nebular gas in front of the central star of the 
compact PN SwSt 1. The absorbing nebular gas exhibits a low degree of ionization, consistent 
with previous studies of this object. We find that the absorbing species in SwSt 1 lie at a 
heliocentric velocity of —40 to —50 km s _1 , in good agreement with Na D measurements. 

2. From the column density ratio of S III fine structure levels, we derive an electron 
density n e = 8800+^ cm 3 , which is at least a factor of three lower than found in previous 
optical and UV emission line studies. 

3. We derive the abundance of iron relative to sulfur, accounting for populations of 
all of the low-energy levels of the relevant ions. We find [Fe/S] = — 0.35±0.12 or —0.50 ± 
0.12, depending on the adopted /-values. The gaseous iron abundance is insensitive to 
uncertainties in n c and T e , and is high in comparison to interstellar medium values, and to 
the Fe abundance derived from optical [Fe III] emission lines in SwSt 1 ([Fe/H]= — 1.64±0.24 
and [Fe/S]= —1.15 ± 0.33). This indicates the presence of marked inhomogeneities in the 
dust-to-gas ratio in different regions of the nebula. 

4. The high gaseous Fe abundance indicates that other refractory elements should not 
be strongly depleted in the sight line. Thus, the elemental Ge abundance lies at the low 
end of the range suggested by SDB. We re-derive [Ge/S] and find a value of 0.72±0.06. 
The large Ge overabundance implies that the progenitor of SwSt 1 experienced s-process 
nucleosynthesis and efficient third dredge-up while on the AGB. 
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5. We do not detect nebular H 2 absorption in the FUSE spectrum. We estimate an 
upper limit to the column density of iV(H 2 ) < 7 x 10 14 cm -2 . This limit is more than 
four orders of magnitude smaller than the beam-averaged column density derived from near- 
infrared vibrationally excited H 2 lines. Apparently, little of this molecular material lies along 
the line of sight to the central star. 

6. The low electron density, non-detection of nebular H 2 , and high gaseous Fe abundance 
derived from the FUSE spectrum of SwSt 1 implies that the gas lying in the direction of 
the central star has markedly different properties from the regions of the nebula sampled 
by emission lines. Since the FUSE spectrum samples only the gas along the line of sight 
toward the central star, this indicates that both the ionized and molecular gas in SwSt 1 
are non-uniformly distributed. However, it is not possible to determine whether the nebular 
inhomogeneities arise from global asymmetries or small-scale clumping of nebular material, 
since SwSt l's angular size is small enough that its internal structure is not resolved in HST 
images. 

7. SwSt 1 exhibits absorption features from excited fine structure levels of neutral oxygen 
at 1040.94 and 1041.69 A. These levels must have been populated by a non-collisional process, 
as indicated by the inverted population structure. Using optical O I line ratios, we find 
that fluorescent pumping by the stellar continuum is the most likely excitation mechanism. 
The fine structure levels of the O I ground term give rise to lines at 63 and 145 jum, whose 
diagnostic value is compromised if the level populations are affected by fluorescent processes, 
rather than being determined purely by collisional excitation. 

The HST data presented in this paper were obtained from the Multimission Archive 
at the Space Telescope Science Institute (MAST). STScI is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. Support 
for MAST for non-HST data is provided by the NASA Office of Space Science via grant NAG 
5-7584 and by other grants and contracts. We are indebted to O. De Marco for making the 
optical spectrum of SwSt 1 available to us, and S. McCandliss for providing H 2 line data. 
We also thank F. Keenan, who supplied us with Fe III level population data, and M. Barlow, 
B. Draine, and S. Federman for helpful suggestions and discussion. T. Kallman provided 
invaluable assistance with the operation of the XSTAR code. This work was supported by 
NASA contracts NAG 5-9239, NAG 5-11597, and NSF grant AST 97-31156. 
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Table 1. SwSt 1 Nebular Component Fits 
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Fe III**** 


1130. 


40 b 


-40.5±1. 


.7 


7 +0 - 4 

' - u -1.8 


1.55±0.25 (13) 


Ge III 


1088. 


.46 


-48.0±0. 


.5 


7.2±1.8 


5.75±0.66 (12) 


I 


1039. 


23 


-58.2±1. 


.5 C 


7.6±2.8 


1.62±0.57 (15) d 


I* 


1040, 


.94, 1304.86 a ' 6 


-43.2±3 


.3 


n-4±g:I 


1.12±0.18 (15) 


I** 


1041. 


.69, 1306.03 a ' 6 


-43.2±3 


.3 




7.94±1.29 (14) 


SII 


1250. 


58 a , 1253.81* 


-43.3±2 


.8 


6.7±1.2 


1.07±0.12 (14) 


S III* 


1015. 


,78 b 


-46.1±0. 


.7 


6.1±0.5 


1.45±0.24 (15) 


S III** 


1021. 


.11, 1021.32 b 


-46.1±0. 


.7 


6.1±0.5 


8.32±1.16 (14) 



Note. - The resulting central heliocentric velocity, Doppler spread 
parameter, and column density are given for the nebular component of 
the indicated species. The fits are displayed as solid lines in Figures 2a 
and 2b. The uncertainties listed are 1-a estimates, and upper limits are 
3-cr. Note that the notation 1.51 (13) = 1.51 x 10 13 . 

a The parameters of these lines are measured from the STIS spectrum. 

b Parameters obtained by simultaneously fitting lines from different lev- 
els of the ion, fixing < v > and b but allowing the column densities of the 
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levels to differ. 

c Highly uncertain < v > due to heavily saturated interstellar compo- 
nents. 

d Uncertain column density due to saturation and/or blending. The 
measured column density is not used in abundance measurements. 
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Table 2. S III Column Density Ratios 



n e 
(cm ) 


-Ro2 a 




4000 

t:W WW 


1 83 


2 87 


5000 


1.42 


2.49 


6000 


1.14 


2.21 


7000 

1 WWW 


96 


2.01 


8000 


0.85 


1.85 


8800 b 


0.77 


1.74 


9000 


0.77 


1.71 


10000 


0.70 


1.63 


15000 


0.51 


1.31 


20000 


0.41 


1.14 


25000 


0.36 


1.04 


31600 c 


0.32 


0.96 


Observed 




1.74±0.38 



Note. - All ratios are computed for an electron tem- 
perature T e = 10, 500 K and the indicated electron density, 
using Shaw & Dufour's (1995) nebular. ionic task in IRAF. 
The adopted transition probabilities are from Nahar (1993), 
while the collision strengths are from Galavis et al. (1995). 

*R 02 = N(S IU)/N(S III**) and 
R 12 = N(S III*)/N(S III**). 

b Adopted n e . Observational (1-cr) errors allow for a range 
of n c from 6400-13600 cm~ 3 . 



c Electron density derived by DM01. Note the poor agree- 
ment with the observed i?i 2 . 
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Table 3. Adopted Sulfur and Iron Level Column Densities 





Index 










Observed 


Adopted 


Ion 


i a 


Level 




Range c 


Rij 


1\ [cm. 


2\d 
) 


S III 







3 p 
-'0 


0.77 


0.53-1.05 




6.41±2.32 


(14) 




1 




3 p 


1.74 


1.36-2.12 


1.74±0.38 


1.45±0.24 


(15) e 




2 




3 p 


1.00 


1.00-1.00 


1.00 


8.32±1.16 


(14) e 


Fe II 





a 


6 n 

-^9/2 


1.00 


1.00-1.00 


1.00 


5.37±0.74 


(13) e 




1 


a 


£7/2 


0.24 


0.16-0.32 


< 0.28 


1.31±0.36 


(13) 




2 


a 


6 n 


0.14 


0.10-0.19 


< 0.29 


7.35±2.44 


(12) 




3 


a 


6 n 


8.8 (-2) 


6.4-12.3 (-2) 


< 0.22 


4.73±1.58 


(12) 




4 


a 


6 n 


4.7 (-2) 


3.5-6.4 (-2) 


< 0.19 


2.54±0.77 


(12) 




5 


a 


4 -^9 /2 


1.02 


0.86-1.08 




5 50±0 97 

/ • / V J — 1 — V J • '7 1 


(13) 




6 


a 


^7/2 


0.15 


0.10-0.22 




7.84±3.24 


(12) 




7 


a 


4 -^5/2 


7.4 (-2) 


4.7-11.7 (-2) 




3.98±1.86 


(12) 




8 
o 


a 


^3/2 


4 Q (-9) 


3 9-7 7 (-9) 

O I.I I ~Zi 1 




2.65±1.22 


(12) 




9 


a 


A E>7/2 


2.7 (-2) 


1.9-4.2 (-2) 




1.46±0.60 


(12) 




10 


a 


4 D 5/2 


9.2 (-3) 


5.6-16.7 (-3) 




4.93±2.99 


(11) 




11 


a 




4.8 (-3) 


2.8-9.3 (-3) 




2.58±1.74 


(11) 




12 


a 


4 D l/2 


2.4 (-3) 


1.4-4.6 (-3) 




1.26±0.86 


(11) 


Fe III 







5 D, 


14.1 


9.9-18.7 


> 3.4 f 


2.82±0.99 


(15) 




1 




5 D 3 


1.64 


1.61-1.64 


1 nc:+0.56f 
i - ud -0.19 


3.28±0.52 


(14) 




2 




5 D 2 


1.00 


1.00-1.00 


1.00 


2.00±0.32 


(14) c 




3 




5 D 1 


0.73 


0.70-0.76 


0.69±0.15 


1.38±0.19 


(14) c 




4 




5 A, 


0.38 


0.32-0.45 


7.8±1.8 (-2) 


1.55±0.25 


(13) c 



a The index i denotes the excitation of the absorbing species, e.g. % = 2 for 
Fe III**. 

b The level population ratio of level % relative to S III**, Fe II, or Fe III**, for an 
electron density of n e = 8800 cm -3 . 

c The range in when n c is varied from 6400-13600 cm~ 3 . 
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d Errors for predicted column densities include the uncertainty in n e , as well as 
those from the measurement of the reference line. 

e The observed column density is adopted for the abundance calculation. 

Saturated line. 
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Table 4. Dependence of Fe and Ge Abundances on Electron Density 



(cm 3 ) 


[Fe/S] a 


[Ge/S] 


[Ge/Fe] a 


6400 


-0.28±0.02 


0.69±0.05 


0.98±0.05 


8800 b 


-0.35±0.02 


0.72±0.05 


1.08±0.05 


13600 


-0.44±0.02 


0.74±0.05 


1.21±0.05 


31600 


-0.58±0.02 


0.76±0.05 


1.39±0.05 


10 5 


-0.68±0.02 


0.77±0.05 


1.53±0.05 



Note. - ogarithmic abundances relative to 
solar (Asplund et al. 2004) are given, for a range 
of electron densities n e . The listed 1-a uncer- 
tainties are not the formal errors, since they do 
not account for uncertainties in n e . 

a Assumes /-values from Morton (2003). If 
NP96 /-values are used for Fe III, [Fe/S] de- 
creases by 0.15 dex, and [Ge/Fe] increases by 
0.16 dex. 

b The adopted electron density and abun- 
dances. 
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Table 5. Optical Iron Emission Line Intensities 





^lab 




io 3 / A / 


Predicted 10 3 / A // H /3 


Predicted 10 3 I x /I w 






Ion 


(A) 


(A) 


Ih/3 


K01 a ; log n c = 4.5 


K92 a ; log n c = 5.0 


(Fe+* 


/H + ) b 


[Fe II] 


8617.0 


8615.9 


0.5±0.1 






2.0 


(-8) 


[Fc III] 


4607.0 C 


4606.6 


1.5±0.1 


1.0 


0.8 


8.0 


(-7) 






^tuo i .u 


i o Q_i_n 7 


19 3 


19 3 


o.o 


( 7) 




4701.5 


4701.0 


5.1±0.2 


5.3 


5.3 


5.2 


(-7) 




4733.9 


4733.3 


2.6±0.2 


2.5 


2.7 


5.6 


(-7) 




4754.7 


4754.2 


2.2±0.2 


2.3 


2.2 


5.3 


(-7) 




4769.4 


4768.9 


1.8±0.2 


1.8 


1.7 


5.3 


(-7) 




4777.7 C 


4777.1 


1.7±0.1 


1.2 


1.3 


7.5 


(-7) 




4881.0 


4880.4 


4.2±0.2 


10.5 


4.4 


2.1 


(-7) 




4930.5 C 


4930.1 


1.1±0.2 


0.7 




8.5 


(-7) 




4987.9 


4986.6 


1.0±0.1 


2.0 


0.8 


2.5 


(-7) 




5084.8 


5084.2 


0.4±0.1 


0.4 


0.3 


5.1 


(-7) 




5270.4 


5269.9 


6.3±0.3 


6.5 


4.0 


5.5 


(-7) 



Note. — Optical Fe line intensities, corrected for reddening using E(B — V)= 0.46 (DM01), 
are given relative to H/3. The optical data were obtained by DM01 from the Anglo-Australian 
Telescope in 1993 May, and were subsequently reduced by those authors. The listed 1-a errors 
account only for statistical uncertainties in the line fits. 

a Predicted [Fe III] emission line intensities (scaled to A4658.1) from K01 = Keenan et al. 
(2001) and K92 = Keenan et al. (1992), for T = 10 4 K and the n c that best agrees with the 
observed intensities. See text for further discussion. 

b The ionic abundance derived from the given line with T e = 10, 500 K and log(n c /cm _3 )=4.6. 
Transition probabilities are from NP96, and collision strengths from Zhang (1996). The adopted 
ionic abundances are given in Table 6. 

c Blended line. 
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Table 6. Iron Abundance From Optical Emission Lines 



Abundance 




or ICF 


Result 


Fe+/H+ 


2.0±0.3 (-8) 


Fe++/H+ 


5.3±0.4 (-7) 


ICF a 


1.16±0.56 


Fe/H 


6.39±3.13 (-7) 


[Fe/H] 


-1.64±0.24 


[Fe/S] b 


-1.15±0.33 



Note. - The adopted ionic and total iron abundances 
are displayed. The errors incorporate uncertainties in n e , 
T e , and the line fits. The ionization correction factor (ICF), 
and hence total Fe abundance, allow for the large uncertain- 
ties in 0/0+ quoted in DM01. 

a The ICF is taken from Rodriguez (2002): Fe/H= 1.1 x 
(0/0+) x (Fe++/H++Fc+/H+). 

b [Fe/S] is derived using the S abundance from DM01. 
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Table 7. Column Density Upper Limits for H2 Lyman Lines 



Line 


A (A) 


/ 




iV(Fit) b 


iV(RMS) c 


8-0 R(2) 


1003.98 


1.662 


(-2) 




< 7.89 


(13) 


8-0 P(2) 


1005.39 


9.911 


(-3) 


< 7.59 (13) 


< 1.40 


(14) 


0-0 Q(2) a 


1010.94 


2.381 


(-2) 




< 4.99 


(13) 


7-0 P(2) 


1016.46 


1.023 


(-2) 




< 9.44 


(13) 


4-0 P(2) 


1053.28 


9.021 


(-3) 




< 1.05 


(14) 


2-0 P(2) 


1081.27 


4.699 


(-3) 








1-0 P(2) 


1096.44 


2.366 


(-3) 


< 4.57 (14) 


< 3.55 


(14) 


7-0 R(3) 


1017.42 


1.838 


'-2) 


< 5.37 (13) 


< 6.07 


(13) 


4-0 R(3) 


1053.98 


1.336 


(-2) 


< 8.32 (13) 


< 5.61 


(13) 


4-0 P(3) 


1056.47 


9.556 


(-3) 


< 6.46 (13) 


< 8.37 


(13) 


3-0 P(3) 


1070.14 


7.537 


(-3) 




< 1.24 


(14) 


2-0 R(3) 


1081.71 


6.363 


(-3) 


< 2.82 (13) 


< 1.43 


(14) 


1-0 P(3) 


1099.79 


2.526 


(-3) 




< 2.90 


(14) 


5-0 R(4) 


1044.54 


1.548 


(-2) 


< 5.50 (13) 


< 4.77 


(13) 


4-0 R(4) 


1057.38 


1.294 


(-2) 


< 3.55 (13) 


< 7.49 


(13) 


3-0 P(4) 


1074.31 


7.744 


(-3) 




< 9.73 


(13) 


2-0 P(4) 


1088.80 


5.151 


(-3) 


< 1.23 (14) 


< 1.98 


(14) 


0-0 Q(5) a 


1017.83 


2.389 


(-2) 




< 4.39 


(13) 


3-0 R(5) 


1075.24 


9.279 


(-3) 




< 9.39 


(13) 



a ine from the Werner band. 

b Upper limit derived from fit within the noise at the neb- 
ular velocity. 

c Upper limit derived from RMS noise in the adjacent con- 
tinuum. 
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Table 8. Estimated Ar tot (H 2 ) Upper Limits for Various Values of T rot 



J 


jVj a 


/j(200 K) b 


/j(500 K) b 


/.I (800 K) b 


/j(1600 K) b 


Mot (200 K) c 


Mot (500 K) c 


Mot (800 K) c 


Mot (1600 K) c 




(cm -2 ) 










(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


2 


< 5.0 (13) 


7.3 (-2) 


1.4 (-1) 


1.3 (-1) 


9.1 (-2) 


6.8 (14) 


3.5 (14) 


3.7 (14) 


5.5 (14) 


3 


< 5.0 (13) 


2.4 (-2) 


2.2 (-1) 


3.0 (-1) 


2.8 (-1) 


2.1 (15) 


2.3 (14) 


1.6 (14) 


1.8 (14) 


4 


< 4.0 (13) 


4.0 (-4) 


2.5 (-2) 


5.6 (-2) 


7.8 (-2) 


1.0 (17) 


1.6 (15) 


7.1 (14) 


5.1 (14) 


5 


< 4.0 (13) 


2.0 (-5) 


1.8 (-2) 


7.3 (-2) 


1.7 (-1) 


2.0 (18) 


2.2 (15) 


5.5 (14) 


2.3 (14) 



a The H 2 column density in rotational level J, from the upper limits in Table 7. 
b The fraction /j of total H 2 molecules in level J, for a given temperature. 
c Total H 2 column density for a given rotational temperature. 
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Table 9. Optical O I Emission Line Ratios 





Observed 


Orion 


Line 


J(A)//(A8446.5) 


J(A)//(A8446.5) a 


I] A4368.2 


1.0±0.2 (-1) 


8.3 (-2) 


I] A5512.8 


4.9±0.5 (-2) 


2.7 (-2) 


I] A5554.9 


4.2±0.8 (-2) 


2.8 (-2) 


I] A5958.5 


8.1±0.8 (-2) 


4.3 (-2) 


I A6046.4 


1.5±0.1 (-1) 


1.0 (-1) 


I A7002.1 


5.6±0.4 (-2) 


9.8 (-2) 


I A7254.4 


4.0±0.5 (-2) b 


1.2 (-1) 


I A7774.2 


3.3±0.5 (-2) 




I A13164 


0.3±0.1 




H I A4861.3 


129 


113 


H I A6562.8 


102 


326 



Note. — Optical line ratios of O I, Ha, and FL3 are 
given relative to O I A8446.5, where /(A8446.5) = 
(2.5 ± 0.1) x 10" 13 erg cm" 2 s" 1 A" 1 . The opti- 
cal data were obtained by DM01 from the Anglo- 
Australian Telescope in 1993 May, and were subse- 
quently reduced by those authors. The l-a errors 
account for statistical uncertainties only. 

a Orion emission line ratios from Esteban et al. 
(2004). These ratios are representative of stellar 
continuum-excited fluorescence. 

b Uncertain - contaminated by telluric feature. 
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i 1 1 r-5 — | 1 1 1 1 1 1 1 1 1 1 1 1 r 



1 120 



11 25 



Fe II I* 




1130 1135 1140 

Wavelength (A) 



1 145 



1 150 



Fig. 1.— A portion of the FUSE LiF 2A spectrum of SwSt 1 from 1120 to 1150 A. The 
data have been rebinned by three pixels, corresponding to a pixel size of 0.02 A. Nebular 
absorption lines (which are much narrower than the broad stellar features) are identified 
above the spectrum; asterisks denote absorption from excited levels. 
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Fig. 2a. — Continuum-normalized absorption profiles are plotted versus heliocentric velocity 
for the identified species. S II AAI250.58, 1253.81 and O I AAI304.86, 1306.03 are from the 
STIS E140M spectrum of SwSt 1. All others are from the FUSE spectrum: Ge III A1088.46 
is from the LiF 2A segment, while the rest are LiF 1A data. Double panels are shown 
for simultaneous fits, and single panels for single line fits. The measured ion is indicated 
above the panel (with asterisks signifying absorption from excited levels), and wavelengths 
are shown within each panel. Model fits to each profile, convolved with the instrumen- 
tal spread function, are depicted as thick solid lines. Thin dashed lines represent fits for 
the individual components before instrumental broadening. Nebular absorption occurs near 
—45 km s _1 . Other species present in the profiles are indicated. See text and Table 1 for 
further information about the displayed fits. 
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Fig. 2b. — Same as Figure 2a, except the displayed continuum-normalized profiles are of 
iron lines from the LiF 2A segment of the FUSE spectrum. 
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0.0 

1.0 



0.0 

1.0 




2-0 P(l) M07B.93 



0.0 

1 



0.0 

1.0 



1-0 R(0) M092.20 




7-0 R(3)_\1017.42 
1.0 -pfljlj ■ 




5-0 R(4) M044.54 

o M^m ,.- /fa>iiJ - jii 



0-0 Q(5) M017.83 




0.0 

1.0 




1-0 P(_l) X1094.05 



0.0 

1 




0-0 0(2)_M010.94 



0.0 

1.0 



0.0 

1.0 



s k^^w^M 1.0 -%*vy^W 



0-0 R(0) M108.13 




4-0 P(3) M056.47 
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3-0 R(5) M075.25 
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Fig. 3. — Continuum-normalized profiles of selected H2 transitions from Table 7 are plotted 
against heliocentric velocity. For comparison, we also display relatively weak (small /-value) 
lines from the J = and 1 levels, for which no iVj upper limits were measured. Two 
transitions from each rotational level are shown, and are arranged from top to bottom in 
order of increasing J. All lines are from the Lyman band, with the exception of the Werner 
lines 0-0 Q(2) A1010.94 and 0-0 Q(5) A1017.83. Absorption from other species is indicated 
within the panels. Fits to the interstellar component (and blends with nearby lines) are 
depicted by thick solid lines. The dotted lines represent nebular absorption profiles with 
column densities equal to the measured upper limits given in Table 7 (if both iV(RMS) and 
iV(Fit) are available, we select the smallest of the two). Dashed lines illustrate the absorption 
profiles expected if the nebular column density of H2 is as high as indicated by the infrared 
emission lines, assuming T rot = 1600 K (for lower T rot , the nebular absorption profiles would 
be even more deeply saturated, since the emitting v = 1, J = 3 level would have a smaller 
fractional population). 



